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Many problems in computational magnetics involve com- 
putation of fields which decay within a skin depth <5, much 
smaller than the sample size d. We discuss here a novel per- 
turbation method which exploits the smallness of e = 8/d and 
the asymptotic behavior of the solution in the exterior and in- 
terior of a sample. To illustrate this procedure we consider the 
computation of the magnetic dipole and quadrupole moments 
of an anisotropic, unconventional, three dimensional super- 
conductor. The method significantly reduces the required nu- 
merical work and can be implemented in different numerical 
algorithms. 



I. INTRODUCTION 

In various problems of electrodynamics, field penetra- 
tion is characterized by a skin depth, 5, much smaller 
than the actual sample size, d. For many purposes, the 
approximation of zero skin depth, or that eddy currents 
flow only as surface sheet currents, is not sufficiently ac- 
curate, as a detailed knowledge of the field penetration or 
accurate values of magnetic multipoles are required. It 
is then important to include the small corrections arising 
from finite skin depth. 

In this paper we examine the inclusion of skin depth 
corrections, focusing on the magnetic response of a high 
temperature superconductor (HTSC). For HTSC's, it is 
knownBia that examining the field penetration yields im- 
portant information about the unconventional electronic 
pairing states in these materials and the still unknown 
nature of high temperature superconductivity. 

Although we focus on this problem, our ideas have 
broad validity and applications. The magnetic response 
of a superconductor is related to that of an ordinary con- 
ductor in an harmonically applied field. The skin effect, 
with 8 <C d, for a quasi-static regime where the frequency 
is restricted by ui -C c/d (c is the speed of light), mapsQ 
onto the problem of a superconductor of the same size 
and shape in a static applied magnetic field. The role of 6 
for a conductor is taken up by A, the effective penetration 
depth of a superconductor (A <C d) . Except for the sim- 
plest geometries where an analytic solution exists for the 
corresponding boundary value problem, obtaining small 
skin depth corrections can be computationally demand- 
ing. These difficulties arise from the nontrivial boundary 
conditions (including open boundary equations at infin- 
ity) and the requirement that the appropriate, generally 
nonlinear differential equations be solved very accurately 



within the narrow region where the skin effects are con- 
tained. Our method offers a way to resolve these difficul- 
ties. 

We examine a superconductor in an applied uniform 
magnetic field, H a . The sample occupies a bounded re- 
gion SlcR 3 and at its boundary, dtt, is surrounded by 
vacuum. For H a smaller than a critical value, a supercon- 
ductor is in the Meissner regime: The magnetic flux is ex- 
pelled from the bulk of the sample. On R 3 \J7 the current 
is j = and it is sufficient to find a magnetic scalar poten- 
tial <&, H = — V$, which satisfies the Laplace equation. 
On f2, the appropriate Maxwell equation is Ampere's law. 
For unconventional pairing states in HTSC's, the London, 
relatione j = j (A) between current-and vector potentials 
A, is nonlinear and nonanalyticElcm Thus it is advan- 
tageous to combine Ampere's law in terms of the vector 
potential with the relation j(A): 
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VxVxA= — j(A). 
c 



(1) 



HTSC's in general have a highly anisotropic structure 
with different penetration depths, Ai, along the various, 
i = a,b,c crystallographic directions. We include this 
penetration depth anisotropy through the anisotropic, 
nonlinear, relation j(A) given in Ref. M. By A we shall 
denote the effective penetration depth (a function of Aj), 
which plays the dominant role in the field decay stud- 
ied. In the special case of an isotropic superconductor 
with a linear relation j(A), all the fields on fl satisfy the 
vector Helmholtz equation V 2 F = F/A 2 , where F can 
be H, j, A. The boundary conditions aroj. — V$ = H a , 
at infinity, while on d£l H is continuousa and there is 
no normal component of current, j n \on = 0. From the 
open boundary condition at infinity combined with the 
the continuity requirement it appears that to obtain the 
finite skin depth corrections, one would have to solve nu- 
merically the appropriate equations in all space. 



II. PERTURBATION METHOD 

To resolve these difficulties, we view the finite skin ef- 
fects, i.e., for finite A in a superconductor, as a small cor- 
rection to the dominant perfect diamagnetic response at 
A = 0. When skin effects are studied, one has to include 
these corrections, which are characterized by the small 
parameter e = X/d <C 1. The boundary value problem in 
the e = limit is relatively simple, one has only to solve 



the Laplace equation for the scalar potential on R 3 \ 17 
with trivial Neumann boundary conditions on <917. We 
assume that an accurate, either analytical or numerical, 
solution on R 3 \ 17 in the e = limit is available. El This 
will be the starting point from which we shall develop 
our perturbation method. The small skin effect is then 
treated as a perturbation from the e = solution. 

To proceed with the perturbation calculation we con- 
sider the auxiliary problem consisting of Eq. (tth on 17, 
the e = solution on R 3 \17, and the boundary conditions 
on 917, j n = and continuity of the tangential compo- 
nent of H (the continuity of H n can not be imposed, it 
vanishes for the external fields in the e = limit). This is 
computationally simple, as it decouples the solutions for 
the regions R 3 L17 and 17. From this auxiliary problem 
we can generated the skin corrections to leading order in 
e, as we shall now see. 

In this paper we consider the magnetic moment of a 
superconductor for an arbitrary direction of H a in a sam- 
ple without a rotational symmetry (this is an extension 
of Ref. 0), and the magnetic quadrupole moment. We 
match the asymptotic behavior of the solution on R 3 \ 17 
and that in 17 by. employing integral identities for mag- 
netic multipolesE At large distances from 17, the multi- 
pole expansion of the fields can be considered and the 
asymptotic behavior is governed primarily by the low- 
est nonvanishing multipole term. There are two different 
ways to obtain the multipole moments: by examining the 
asymptotic behavior on R 3 \ 17, and from the fields com- 
puted on 17. By matching the asymptotic behavior we 
mean that the exact solution on R 3 \17 is formally written 
in terms of the unknown multipole moments which must 
agree with those computed from the fields on 17. We for- 
mulate integral identities to compute the magnetic mul- 
tipoles from fields on 17 such that we can identify terms 
in these expressions which are of different orders in e. 

The magnetic momenta is 



m = 



1 /dl7r"xj(r"), 
lc Jn 



(2) 



where r" is the position vector for a point in 17 and j is 
found from Eq. (El). The components of m can be writ- 
ten for e C 1 in the form m, = TOoj(1 — cti e + 0(e 2 )), 
i = x,y, z, where m-oi, denoting the i — th component of 
the magnetic moment in the limit e = 0, represents a per- 
fect diamagnetic response. For an ellipsoid it is given by 
a demagnetization factor. The oti describe small correc- 
tions to perfect diamagnetism due to current penetration. 
For a direction i, where m oi — (it could vanish from 
symmetry arguments, for a particular direction of H a ), 
one can still have qjj ^ 0, because of the anisotropic and 
nonlinear relation j (A) . The effects of nonlinearity j (A) , 
absent for e = 0, are typically small and can be thought 
of as field dependent corrections to aj, linear in H a . 

To distinguish terms in Eq. (0) of various orders in e, 
we use Ampere's law, identities from vector calculus, and 
Gauss' theorem to obtainll 



m = Z~ I dS l n ( r " • H) + n x (r" x H)] (3) 

87r JdQ 
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gJ17 H = mi + ni2, 



where r" is the position vector for a point on 917 and 
n is the unit normal pointing outwards. The terms mi 
and rri2 are of different order in e and the latter is small, 
i.e. of O(emo). This can be seen from the expression for 
ni2. Since H is confined to a "skin" layer of thickness A, 
the integral over the whole volume of 17 is effectively only 
an integration over the region ~ A away from its surface. 
Thus rri2 vanishes in the zero penetration limit (e = 0) 
and m(e = 0) = m = ni!. In order to obtain m to 
O(emo) it is sufficient to compute ni2 to leading (zeroth) 
order. The term 1112 explicitly scales with e and any 
first order corrections for the fields needed to compute it 
would only produce contributions of order O(e 2 mo). 

A similar integral identity can h£ derived for the mag- 
netic quadrupole moment, defineaj as a symmetric trace- 
less tensor with components 



Q ij = ijf ( fil[r"(r"xj) + r"[(r"xj)] 
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(4) 



Using the previously introduced notation for terms of 
different order in e, Qij = Quj + Q2ij we can derive, 
employing integration by parts and standard identities: 



Qui — 7— 



8tt Jan 
\r"H. 



fr'^Kn-r")], 



(5a) 



Qa« = ^- I dSl [3 [r'/H, + rj'fli] - 2%(r" • H)], (5b) 
S7r Jn 

where <5y is the Kronecker symbol. If by Qoij /Owe de- 
note a particular component of the magnetic quadrupole 
tensor in the e = limit, then as in the case of the mag- 
netic moment, we conclude that Quj is of O(Qoij) while 
Q 2 ij is of O(eQ 0ij ). 

Using Eqs. (0) and (0) we can match the asymptotic 
behavior of solutions in regions R 3 \ 17 and 17. It is 
then possible to perturbatively obtain physical quanti- 
ties to leading order in e, using only the fields on 17 
computed from the auxiliary problem. The fields and 
quantities evaluated from this problem are denoted by 
an overbar notation. We consider first the magnetic mo- 
ment. On R 3 \ 17, the scalar potential can be written as 
$ = $ a + $ r , where $ a is the potential due to the ap- 
plied field and satisfies the open boundary condition at 
infinity, — V$ a — > H a , and $ r describes the presence of 
the superconductor. Since m2, as we have shown, explic- 
itly scales with e, it can be accurately computed to first 
order in e by obtaining its leading contribution. This 



is achieved by using the fields on fi from the auxiliary 
problem, i.e., by writing $ = <f> a + $ r . The task of de- 
termining m to O(emo) is therefore reduced to that of 
correctly including the contribution of mi to first order 
in e. The exact solution for H is continuous on <9f2. To 
calculate mu, the j-th component of mi, we can use the 
external fields obtained from $. The part of <& r which 
has a dipole character is characterized by the unknown 
vector m, the correct value of the magnetic moment. The 
remaining part of $ r has, different symmetry properties 
and does not contributed to TOj. Contributions to the 
i-th component of mi, i = X,y, z, (we take m$i ^ 0) can 
be written as 



by incorporating our perturbation procedure in the ap- 
propriate numerical algorithm. The integral identities 
for the magnetic multipoles are valid in the quasi-static 
regime and not restricted to the field of superconductiv- 
ity. 

We have addressed here the computation of small 
skin effects in a superconductor, using the perturbation 
method and matching the asymptotic behavior of a so- 
lution. We have shown how to accurately compute the 
nonlinear magnetic response of an anisotropic supercon- 
ductor, by simplifying the boundary conditions and re- 
ducing the size of the computational domain. 



mij(m) = mii($ a ) + m u (^ r ), 



(6) 



ACKNOWLEDGMENTS 



where mu(<fr a ), for $ a , which is known, can be simply 
calculated from Eq. (Q). We define py by mu($a) = 
J2jPij m oj7 hj — x,y,x. The constants py, which de- 
pend on the shape of f2, can now be determined by solving 
for pij using the known values mii(& a ) and m Qi corre- 
sponding to H a applied along three independent direc- 
tions. In the limit e = 0, vti2i — and from Eq. (|6j) 
follows the identity 



m i 



mu 



E 



Pij m 0j + ^2 (% _ Pv ) m oj > ( 7 ) 



where mij($ r ) = X);(% ~ Pij) m 0j- F° r c / 0, when 
the solution for <& and H is given in terms of a multi- 
pole expansion with unknown coefficients, mu(^ a ) re- 
mains the same. The terms in $ r which contribute to 
the magnetic moment m, will now have coefficients pro- 
portional to the correct unknown value of m, slightly 
changed from the e = case. We can therefore write 
mii(& r ) = J2j($ij ~Pij) m j an d the correct value for mu 



satisfies mu = J2-j[ r 



P t j( r 



moij)]. Employing 



the fact that m-n and m-n agree to O(emoi) we can solve 
for m l from m, = Y^j[^ij m ij ~ Pi 3 { m ij ~ m 0ij)] + m 2l , 
with the solution for rrii correct to O(e), 



m 0l 



J2Pij lm 2j 



m 0i 



^P i3 l m 2 j. (8) 



Therefore the magnetic moment can be computed by only 
determining the lowest order contribution to ni2. 

Following an analogous procedure we can obtain a so- 
lution for the components of the quadrupole tensor Qy 
accurate to first order in e. The resulting expression is 
similar to Eq. (||) with the constants py replaced by the 
appropriate fourth rank tensor. 

As a simpler example, it is instructive to consider an 
isotropic, linear superconducting sphere Ja an applied 
uniform magnetic field. Eq. (g) reducesEJ to the ana- 
lytical resultE m = mo(l — 3e), where mo = —H a a 3 /2, 
e = A/a and a is the sphere radius. 

The ideas presented here can be used in problems in 
computational magnetics involving a small skin depth, 
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